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Resonances of cylindrical symmetric microwave cavities are analyzed in 7?.-matrix theory which 
transforms the input channel conditions to the output channels. Single and interfering double 
resonances are studied and compared with experimental results, obtained with superconducting 
microwave cavities. Because of the equivalence of the two-dimensional Helmholtz and the stationary 
Schrodinger equations, the results present insight into the resonance structure of regular and chaotic 
quantum billiards. 
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I. INTRODUCTION 

Nuclear reactions at low energies are characterized by narrow resonances which, in the case of well separated 
resonances, show a Breit-Wigner shape. If the separation is small, however, interference occurs which can lead to 
complicated non-homogeneous resonance shapes. This has been observed already in the early days of the study 
of slow neutron resonances. These data stimulated the theoretical development of resonance theories, as e.g. the 
stationary 1Z- matrix theory of Wigner and Eisenbud Q] , or the equivalent one for running waves of Kapur and Peierls 
■ The idea behind this work is a transformation from the input to the output channels without taking into account 
explicitely the complex dynamics of the strongly interacting system, which enters only through the parameters of the 
-h ■ transformation matrix. These parameters have to be obtained by fitting the resulting multi-channel cross-section to 
experimental data. The essential benefit of the Wigner and Eisenbud formulation is the transformation to standing 
waves which results in a spatially stationary problem with real boundary conditions rather than the complex boundary 
£NJ ' conditions of running waves in each channel. Nevertheless, this procedure gets rather complicated when more than 
one resonance contributes significantly, and empirical analysis can be hindered considerably by insufficient statistics 
and large background contributions. A complete representation of the theoretical as well as the practical aspects 
of 7?.-matrix theory in nuclear physics which constituted the basis for the widespread later work in this field is the 
review article of Lane and Thomas |3( ■ Although in a different context than in the present paper, 7?.-matrix theory 
has recently also been applied to multiphoton processes 0, , in the analysis of Rydberg spectra 0, , and to chaotic 
£T) . ionization of Rydberg states 0, @ . 

A sensitive test of single- and multi-level resonance theories can be obtained in the electromagnetic analogue of a 
quantum mechanical multi-channel system in the form of a flat superconducting cavity, fed with microwave radiation 
[loL ITU . It has been shown [Tjj that all resonances of such a system can be resolved and analyzed with respect to 
their positions and widths. This allows a clean verification of level statistics for regular and chaotic cavities, resp., 
depending on the geometrical shape of the flat resonator. Superconducting cavities, with their high quality factors, 
can, however, also be analyzed with respect to line shapes of isolated and overlapping resonances which are almost 
undisturbed by additional broadening from wall losses leading to strongly overlapping resonances in the spectra, as 
in normal conducting resonators. This is typically due in the high frequency part of the spectra, relevant especially 
for chaotic systems [[]]]. There, the 'true' resonances, i.e. not additionally broadened, still have widths smaller than 
the mean level distance which allows a theoretical approach on the basis of groups of isolated resonances. First 
analyses of resonances with the Breit-Wigner formula have been performed in j 1 12L [l3j . We amend this work in two 
aspects. First, we start from the electromagnetic field conditions of microwave cavities, rather than assuming the 
Breit-Wigner formula as valid also there, and second, we use the 7?.-matrix formalism, rather than starting directly 
with the iS-matrix, since the parameters of the 7?.-matrix allow direct modeling of resonance shapes and comparison 
of resonances in different channels. 

In this paper, we first derive the resonance theory of a two-dimensional microwave cavity, following the idea of 
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Wigner and Eisenbud |lj to parametrize the relation between input and output in the form of a multi-level reactance 
matrix. From this, the S'-matrix, giving the power transfer from input to output, can be constructed by a non-linear 
transformation which leads, in the multi-level case, to rather involved interference structures. 

The approach is tested with experimentally determined single and double resonances of a superconducting two- 
dimensional billiard of threefold (C3) symmetry ^[^jj. We also analyze single resonances of a three-dimensional 
stadium billiard 0|. The tests are performed in the low as well as in the high frequency parts of the spectra, in 
order to cover, on one hand, strong resonances which extend over several orders of magnitude, and on the other hand 
the relevant frequency part for analyzing the dynamics of the system. At high frequencies level interference occurs, 
leading to characteristic non-uniform line shapes. Comparing our one-level formula with the Breit- Wigner formula of 
nuclear physics shows that for narrow and isolated resonances the difference is neglibly small. 



II. DERIVATION OF THE CAVITY RESONANCE FORMULA 



A. The cavity and antenna wave functions 



We treat a cylindrical microwave cavity with N antennas 1, ... a, ... /3, ... N connected perpendicular to it. One 
antenna, which we always denote by a, is feeding the cavity while radiation is extracted through the others (multi- 
channel situation, see Fig. Q). The electromagnetic fields in the cavitiy can be derived from scalar fields 0]> and it 
suffices to treat these. The complete scalar wave function with fixed frequency ui can be written as 

S(f, t) = *(r) exp(-iut) ; r = {x,y,z} (1) 

with the spatial part ^(r) obeying the wave equation 

(A+^)*(f)=0. (2) 

For a cavity with arbitrary cross-sectional shape, but with end-plates perpendicular to the cylinder (cf. Fig. [JJ, the 
spatial part 'i'(r) can be separated into longitudinal and transverse components, according to 'i'(r) = ip(z)x(r't) with 
z, the longitudinal and f t = x ■ e± + y ■ e*2, the transverse coordinates. The z-dependent part consists of standing 
waves, A sin kz + B cos kz, and the boundary conditions, E t = at z — and z = d for transverse magnetic (TM) 
fields lead with k = pir/d to the form 

J)7T Z 

E z =ip T M,p( z )x(x,y) = cos(—)x(x,y), p = 0,1,2,... (3) 

while for transverse electric (TE) fields the vanishing of H z at the z-boundaries requires 

P7T z 

H z = ip TE ^ p (z)x(x,y)=sm(—)x(x,y), p = 1,2,3,.... (4) 

For flat cavities, i.e. d much smaller than the wavelength of the resonating microwaves in the cavity, TM modes 
are the only possible ones since the p = modes alone can exist. This is the interesting case for a comparison 
with quantum mechanics [Tl| . However, the following derivation is more general and holds for any value of p, thus 
describing arbitrary cavities with cylindrical geometry. The transverse part x(r) is a solution of the two-dimensional 
wave equation 

(A t + f) Xp (r t ) = 0; A t = A - ^ ; = ^ f-ff (5) 
with the boundary condition on the cylindrical surface S (cf. Fig. ^) 

X P \s = 0- (6) 

Equation Q together with Eq.©, defines an eigenvalue problem for the two-dimensional transverse fields with 
eigenfunctions Xpr(^) , T = 1, 2, 3, . . .. For p = 0, these fields are the direct analog of the stationary wave functions of 
a two-dimensional quantum billiard which makes the study of the flat electromagntic cavity interesting for illuminating 
the implications of quantum chaos. The TM fields are given by 

pir c 2 ,PTtz - 

Q IWC pTTZ 

Ht = — J- C0S(— — jVtXpr- 
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The whole system can geometrically be divided into an internal and an external part where the internal part 
consists of the cavity proper, while the external parts are the N antenna wave guides. The concept of the 7?.-matrix 
approach consits of obtaining a transformation from the feeding antenna to the output antennas by using the boundary 
conditions at the connections between internal and external regions. 

The longitudinal wave in an arbitrary antenna fi out of the N antennas in general consists of incoming and outgoing 
parts 

ipniZfi) = A n exp(-i/c M z AI ) + A'^ exp(ik^z^) . 
The transverse part x^if^t) obeys also the reduced wave equation 

,.,2 Q2 

For the group velocity of the running wave in antenna /i one obtains 



(A fit + ^)xM = 0; = A — j w£ = w a - (7) 



10J 



dk.j 



(8) 



We now specify to the case of N open antennas where the input into antenna a produces outgoing fluxes in a 
(direct reflection) as well as in all other antennas [3 ^ a. For this purpose we introduce local antenna coordinates 
for an arbitrary antenna (i according to 

X H — x ~ x fj, i Vy, = V ~ Vn i z fi = z — z n 

and {x^q, j/ M 0j z ^o\ denotes the footpoint of antenna fx. For the geometry of Fig. ^ we have z^q — zq — d for all N 
antennas. The spatial wave function, \&(f), depends on which of the N antennas has been chosen as the input antenna, 
and thus contains incoming waves. For fixed frequency u) there are N such distinct wave functions which form a basis 
for the scalar cavity fields. To distinguish them they have to be marked by the index of the feeding antenna. With 
this notation the basis set consists of the N functions $ M (r) ; fj, = 1 • • • N. With our convention, which alwas denotes 
the input channel with a, the spatial wave ^ a (r) approaches in antenna (3 of the external region the form 

* Q (f— ► rp) = il> a p{zp )xpl?pt) inantenna/3. (9) 

Introducing the unitary S'-matrix which transforms input into output we can write the longitudinal parts in the 
form 

ipapizp) = Sap exp(ikpZp)Kp; [3^ a 
ip aa {z a ) = [exp(-ik a z a ) - S aa exp(ik a z a )]K a 

= (ipin + Vw) K a (10) 

where Kp, K a are normalizing constants. For later purposes it is convenient to choose normalizations such that the 
ingoing flux per unit time is normalized to one, and the outgoing fluxes per unit time are given by the absolute squares 
of the S'-matrix elements. This leads to 

= Kl V4V>m ' J X*aXad(T a ■ V ga = 1 ; J X* a Xad<J a = f a 

a 

K a = .. = l(-^)Va (11) 

P a r = K^f^p 1 \S a pf ; Kp = I ^ f2 = 

p out 

P a ^t = M ^ _ = | 5a/? |2 (12) 

a 

where a is the input antenna, (3 ^ a are the N — 1 output antennas, and da a is the transverse surface element in 
antenna a. 
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B. Transformation to standing waves 

Analog to the formalism of Wigner and Eisenbud Q we perform a basis transformation to standing spatial waves 
$ a (r), where the index a again indicates the input antenna, and the whole set {^p , (3 — 1 ■ ■ ■ N} is also complete for 
fixed uj. This has the benefit to render the boundary conditions real. Because of the completeness of both sets we 
can expand the function $ Q into the set p\ 

®a = ■ ( 13 ) 

P 

The new base function <fr a (r) assumes, in analogy to Eq.(9), in the external region the form 

$ a (f-> rp) = <f> a p(zp) xp(rpt) inantenna/3 (14) 

with 

4>ap{zp) = COS (fc,3 2/5)^/3 P ^ a 

<t> aa (z a ) = [K 1 sin(k a z a ) + K aa cos(k a z a )]Z a , (15) 

where 1Z is the derivative (Wigner and Eisenbud 0), or reactance (as in most work on electromagnetic applications), 
matrix, and a factor k~ x has been put in for later convenience. From Eq. l|15|l follow the derivatives at the antenna 
footpoints 

(^)|o=0; A£L)| =Z o . (16) 

ciZ/s az a 

We can visualize this, see Fig- El by placing virtual mirrors (Mi) into the antenna wave guides such that the system 
becomes stationary with standing waves in each antenna. The mirror is located in all exit antennas a quarter 
wavelength away from the antenna footpoint, while in the feeding antenna it has to be placed in such a way that the 
frequency of the system assumes the freely choosable value u. 

Now, in the external region, we can identify the coefficients of the independent functions exp(i/c At z /J ) and exp(— ik^Zp) 
on both sides of Ea. (|13fl with the explicit forms of the wave functions in the external region, Eqs. i|10fl and (|14[) . to 
arrive after a straightforward calculation at the connection between 1Z- and 5-matrix which in matrix notation reads 



with the diagonal matrices 



and with 



S = p(l-iBKB)- 1 (l + iBKB)p (17) 

B af 3 = 5 a 0kl/ 2 ; p a/ 3 = S af3 exp(-i k a z a0 ) , (18) 



Z, = k^K» = kf)V\ (19) 



Equation l|17|) defines the nonlinear relation between the S- matrix and the reactance matrix B1ZB. 



C. Derivation of the reactance matrix 



We now determine the reactance matrix from the relation between cavity eigenfunctions and antenna wave functions. 
This procedure is the essential part in the Wigner-Eisenbud approach, since it relates the S'-matrix, Ea. (|17|) . to the 
eigenstates of the cavity, without explicitly detailing its properties. We start with expanding the standing wave 
function $ a inside the cavity into the complete orthonormalized basis of cavity eigenfunctions which are the cavity 
eigenstates when all antennas are removed (closed internal region) 

s 

and from Eqs. 10 El EJ) the set \l/ s is defined (for transveres magnetic waves) by 



* s (r) = *Ptm, p ( z ) X P r(x, y); s = {p,r} ; p = 0, 1, 2, . . . ; r = 1, 2, 3, 



(21) 
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where p enumerates the longitudinal TM waves, and r the transverse parts. The orthonormalization implies (though 
in our case of standing waves and time reversal symmetry all functions are real, we denote for generality the adjoint 
functions as in a complex function space) 

J ^ s ,(r)V s (r)d 3 r = 8 S , S , 

and the boundary conditions are (see Eqs. 

Xr P (x,y) \s= 0, d ^™' V \z=o,d = 0- 
The eigenfunctions f s obey the eigenvalue equation 

uj 2 

(A+^-* s =0 (22) 

with eigenfrequencies u) s , while the antenna wave functions have a continuous frequency spectrum 

j 2 

„2 



(A + ^)$ Q =0. (2:-!i 



From Eqs. J20I we obtain 

^{co 2 - u 2 )D as = J {$* A* s - #*A$ a }d 3 r . (24) 



Applying Green's theorem 



, . , . / , dv du 

(u/Xv — vAuja r = / (it— v——)a<j 

On on 



we can convert the right side of Eq. I)24|l into a surface integral which has contributions from the surface S and from 
the horizontal plates z = and z = d: 



2 = 



z—d 



The first integral vanishes, because on S both, <£> Q and \l/ s , are zero according to the boundary conditions. The second 
integral also vanishes, since either s /dn or d<& a /dn are zero on the (z = 0)-plane. The only non-vanishing contri- 
butions come from the third integral at each antenna footpoint(/x = a . . . N). Because of the boundary conditions, 
Ea. (|16|l . we obtain 



l(c 2 -^)^ QS = - W \ oX ,(x„y,)d X dy 



= --(f) 1/2 / M fe/. P (rf) 
c fa 

fsa = / Xi( X a,ya) ' Xa(x a : Va)dx a dy a , (26) 



since only the antenna a — a contributes to the sum. From this we obtain for the expansion of Eq. (121)11 

1 1,0. 
UJ 2 - to 2 c f, 
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Contracting each side of Ea. (|27(l for an arbitrary antenna j3 with the transverse antenna wave x% we obtain from 
Eq.JTSJ, Eq.d23) and with f* p = J XpfapiVp) ' Xs(xf3,yp)dxpdyp 







J Xp^apdxpdyp = K a p^(j-) 1/2 J XpXp dxp dyp 



c 

TJTT* -ST ) V2 /-^ ^TMAd)} 2 ■ (28) 



With 



we finally get 



lsp = -^fifan IpTALs(d) (29) 

n a0 = £ ■ (30) 

s a 

which is determined by the values of the cavity wave functions at the antenna footpoints. If one uses the pro- 
portionality of energy and frequency in quantum mechanics {E = fko) the difference of the reactance matrix 1Z in 
the electromagnetic and in the quantum case [l| is the resonance denominator. Here it contains the squares of the 
frequencies while in quantum mechanics the frequencies enter linearly. 

The properties of the S'-matrix for time reversal invariant systems, unitarity and reciprocity, which lead to relations 
between the S'-matrix elements which are not linearly independent, enter in a simple way into the 1Z- matrix: 1Z is a 
real and symmetric matrix which (for N antennas) has N(N + l)/2 independent real parameters. Since for chaotic 
billiards one eventually wants to study systems which are not time reversal invariant, we have presented the derivation 
for the 7?.-matrix more generally, without using the reality condition. 



III. THE EXPERIMENTAL RESONANCE SPECTRUM 



As pointed out in Sec. II. flat microwave cavities with cylindrical symmetry are a powerful tool for experimental 
studies of two-dimensional quantum billiards. While early experiments were carried out with normal conducting 
devices at room temperature [la . Il9i | , only superconducting resonators allow measuring complete spectra in a large 
frequency range with a high signal-to- noise ratio 0, 0] . 

In [14j a chaotic billiard possessing threefold (C3) symmetry was studied experimentally. Its spectrum shows single 
as well as double resonances (Fig.|3 upper part, and Figs.0]andEl and thus provides a basis to test resonance formulae 
for non-degenerate and nearly degenerate modes. The resonator used in the experiment is made from lead-plated 
copper and becomes superconducting below a critical temperatur of T c = 7.2 K. The height of the cavity is about 
5.8 mm and therefore two-dimensional wave propagation is ensured up to a critical frequency of about 25.5 GHz. The 
experiment was carried out inside a liquid helium cryostat under stable conditions, where at a temperature of 4.2 K 
the power transmitted from one antenna which emits microwaves to a second antenna which picks up microwaves, was 
measured for frequencies from 45 MHz to 25 GHz in steps of 10 kHz. The antennas consist of a wire that penetrates 
the cavity through small holes in its lid by less than 0.5 mm. Due to the low rf resistance the superconducting 
resonator possesses a quality factor Q of the order of 10 5 and hence the lifted degeneracies of the double resonances 
can be resolved. Even higher quality factors (Q 10 7 ) can be achieved when using electron-beam welded niobium 
resonators (T c = 9.2 K) or three-dimensional microwave billiards (see 01)- The latter ones, however, do not show 
an analogy between Helmholtz and Schrodinger equation or cylindrical symmetry, but wave chaotic phenomena of 
the vectorial Helmholtz equation. Such a billiard was studied in |16( and we analyzed here some of its resonances for 
comparison (Fig. [31 lower part, and Fig.^J). 

For our analyses we chose resonances in the low and high frequency parts of the spectra (Figs. |3J H| [5J andEJ). We 
also measured one single resonance of the C3 billiard at room temperature, i.e. with the normal conducting resonator. 
This clearly shows a broader resonance line and a reduction of the power transmitted by several orders of magnitude 
(Fig. 0} . A more detailed description of experiments with superconducting microwave billiards can be found in 111. It 
should be emphasized, however, that, though the measurement of a whole spectrum typically takes one or more days, 
the data points in the vicinity of an individual resonance are taken within a few seconds under identical circumstances. 
This ensures the high precision of our experimental resonance lines. 
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IV. ANALYSIS OF SINGLE AND DOUBLE RESONANCES 



For the analysis of experimentally determined resonances we augment the 7?.-matrix, Ea. (|30() . in two aspects. First, 
we allow for the fact that the superconducting cavity is not a perfect conductor for electromagnetic waves. This 
results in an overall damping of the fields which are no longer strictly harmonic, as assumed in Eq.Q. We account 
for this by adding to the frequency a small imaginary part, u s — > oj s — i\ s , which results in an additional broadening 
of the resonances (cf. [l^, Chap. 8. 8, where cavity losses because of finite wall conductivity are treated explicitely). 

The introduction of a damping factor A produces a redundance problem for the extraction of the resonance pa- 
rameters from the fitted spectra. Since the 7?.-matrix contains for two antennas and s overlapping resonances 3s 
independent parameters (the reduced widths, r y s i,s2, and resonance frequencies, cj s ), we introduce s parameters more, 
the damping factors X s . This ovcrdetermines the analysis since one can shift in the parameter space between the par- 
tial widths, T s i = 7jj, determining the resonance strengths and widths, and the damping factors A s , without changing 
the fit. We have carefully analyzed this problem analytically and numerically. The result is that the uncertainty 
introduced thereby does not influence the leading order of the partial widths and damping factors. 

For the numerical analysis we set c = 1 which then measures the wave numbers k^ and the 7^ in s _1 . 



A. Single resonances 



For single resonances the power transfer between one input and one output antenna, P\i =\ S12 | 2 , is determined 
by Ea. 1)12(1 and Eg. (|17j) which connects the reactance matrix with the ^-matrix. For a single resonance, viz. only one 
term in Ea. l|3U[) . this leads to a resonance formula which resembles in its structure the Lorentz line (for the relation 
between this and the Breit-Wigner form, cf. Sec. V) 

P (uj) = 4 7i 2 fci7 2 2 fc 2 , , 

121 ' - + A 2 ) 2 + ( 7l 2 h + 7 2 k 2 + 2Auv) 2 [ ' 

where in the one level evaluation of Ea. (|30fl the index s = 1 has been obmitted. The partial widths, Ti, and the 
additional damping factor, e, are (the reason for this choice is presented in Sec. V) 

r i = 7 ?; e = 2A (32) 

which results in the total width 



Ttot = Ti + T 2 + e , (33) 

and the lifetime of the resonance is given with the complex resonance frequency as r = A -1 . 

The evaluation of the parameters in Eq. I|31|) for a given cavity resonance is performed with the software package 
Mathematica 4.1, 'Nonlinear Fit'. For the antenna wave numbers ki we use the common value fej = u) (i = 1, 2) which 
is the dominant mode for a coaxial cable |l7j . 

Figure |3] shows two isolated cavity resonances (dotted points) and their analysis with the theoretical relation, 
Eq. (|31|l . (continuous line). The fit to the experimental points over a total interval of 4 to 7 orders of magnitude is 
nearly perfect, and this shows the reliability of the derived resonance formula. The resonance parameters are listed 
in the figure caption. 

Figure 0] depicts the comparison of an isolated resonance in a 'warm', normally conducting cavity with the same 
resonance in the superconducting cavity. The drastic increase of the width is due to the additional wall damping. 
In the analysis this reflects itself in the vastly increased damping factor, e. There is also a shift of the resonance 
frequency, because the superconducting cavity at T=4.2 K is contracted in size, and thus its eigenfrequencies are 
higher (cf. [2(j). For the resonance parameters, see the figure caption. 



B. Double resonances 



In this case the 7?.-matrix is given by (see Eq. I|3U|) ) 

„ 711712 . 721722 r „.v 

y ^ 12 - —2 — -72 + —2 — -p. ■ 

The resulting expression for the S'-matrix, Ea. (|17|) . is extremely involved, containing a lot of nonlinaer interference 
terms, and it can not be simplified since no small parameter for an expansion can be defined. So, we present the 
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resulting formula for P12 only in the Appendix. It is obtained by inserting Ea. l|34|) into Ea. (|17fl . and this into 
Eci. (|12[) . Figure [SI shows two interfering double resonances and their analysis in 7?.-matrix theory. The two resonances 
exhibit very different interference stuctures which both are well described by the 7?.-matrix approach. The resonance 
parameters for both resonances are again listed in the figure caption, where the individual parameters for the two 
resonances have restricted significance because of the strong interference between both. 

In Fig. El we present a section of the high frequency part of the spectrum of the 3D-cavity (from Ref. |l6j |. more 
than 5000 resonances below this section!). It can be seen that the resonances are still well separated, leading to single 
and at most double resonances. The insert presents a fit to the double resonance around 13.033 GHz and shows the 
applicability of our analysis also in this part of the spectrum of a cavity which is not of a cylindrical symmetric shape. 

In all the analyses shown, the quoted fit accuracy \ 2 refers to the performed logarithmic fit of the logarithmic data. 

V. COMPARISON OF THE CAVITY RESONANCE WITH THE BREIT-WIGNER FORMULA 

For the comparison of the two 1-level formulas we start from Eq.Jjnj, i.e. 

4 7l kll2 k 2 



(co 2 - uo 2 r + A 2 ) 2 + ( 7l 2 fci + 7 2 fc 2 + 2Xcj r ) 2 
With ki—k^—ijJ and oj 2 ^ := ui 2 — A 2 one obtains 



p / v = 4 7i 2 7 2 cj 2 

' I2 "- {ijJ- uj eff ) 2 (uj + u! e ff) 2 + (2Aw r + 7 2 w + 7|w) 2 ' 



For narrow resonances r,A<w r the response is non- vanishing only around u) « to r . Then we can set 

to = Lo e fj + 61LU u> = u> r + £>2W, and , <C 1 . 

Up to first order in Siu/uj and ^w/w, Eq. I|35|) . takes the form 



P12H 



7i 2 7l 



(u c e// ) 2 (l - 4^) 2 + I[2A(1 - + 7? + 7 2 2 ] 2 



Neglecting the small corrections Siuj/cj, ^w/w, and with the partial and total widths of Eqs. E3J , Pi2{u) assumes 
the Breit-Wigner form 

P ^ H « 7 r ! 2 r ', lr yi ■ (36) 

To give an example, for the first resonance of Fig. we have T tot /2ir — 294.9 kHz, f r — u r /2Tr — 4.63 GHz and 
thus with Su> « T to t, ^tot/^r = 6.5 ■ 10~ 5 . Consequently, the Breit-Wigner formula is an excellent approximation for 
narrow single resonances in the electromagnetic case. 

VI. CONCLUSION 

In this paper we study the resonance spectrum of a flat electromagnetic resonance cavity which is fed through an 
input antenna and analyzed by one or several output antennas. The TM eigenstates of the cavity obey the same 
Helmholtz eigenvalue equation as the two-dimensional stationary quantum problem in the same geometry (quantum 
billiard). Thus the study of cavity resonances provides, especially in the almost lossless superconducting case, an 
excellent tool to study experimentally the spectral properties of regular and chaotic quantum billiards. It is therefore 
of great importance for the analysis of the elctromagnetic case to derive a reliable resonance formula which not only 
reproduces the resonance shapes, but from which one can also deduce the relevant resonance parameters. 

We derive the cavity resonance formula in close analogy to the Wigner-Eisenbud 7?.-matrix theory, which is a corner 
stone in nuclear reaction theory, and has been used in various fields, especially recently in atomic physics in the 
analysis of Rydberg atoms and multiphoton ionization. The goal of this approach is to express the S'-matrix in terms 
of properties of the wave function at the channel entrancies which in our case are the antenna footpoints at the cavity. 
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We succeed to express the cavity- antenna couplings in terms of the cavity eigenfunctions at the antenna entrances. As 
in the nuclear case, the resonance parameters can then be extracted by fitting measured resonance shapes. Differences 
to the nuclear case come from the different dispersion laws of the electromagnetic and quantum cases, and from the 
geometrical properties of TM resonances in the flat cylindrical cavity. For single, narrow resonances, however, the 
usual Breit-Wigner formula is an excellent approximation to the electromagnetic Lorentz lineshape. 

We apply our multi-antenna result to measured well resolved single and double resonances of a superconducting 
cavity with one input and one output antenna. The fits are almost perfect over many orders of magnitude of the 
transmitted power. Rather involved interference patterns of double resonances (cf. Fig. [SJ are well accounted for by 
the two-level formula, even in the high frequency parts of the spectra. 
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APPENDIX: THE FIT FORMULA FOR TWO ANTENNAS IN THE TWO-LEVEL CASE 

The power transfer P\ 2 from the input to the output antenna is determined according to Ea. (|12j) by the absolute 
square of the ^-matrix element 

pout 

Px2 = 9fr =1 s 12 1 2 . (A.i) 

The S'-Matrix is given according to Eq.(17) as 

S = p{l - iBKB)- 1 ^ + iBKB)p . (A.2) 

Since the diagonal phase factor p does not contribute to the absolute square we omit it in the following. Then we 
obtain with BUB = M 

s (u) = 2iM 12 M 

12V ; 1 - iMu(u)) + Mi 2 (w)M 2 i(w) - !M 22 (w) - M 11 {uj)M 2 i{uj) ' 

In order to simplify the evaluation we absorb the diagonal B-matrix in the partial widths 75, and with ki—io (cf. 
Sec. IV) we set 



7B,yH = Vw7tf5 (*,i) = (l,2). (A.4) 

Now we obtain from Eq.(34) for BUB 

m «.j - ~~a — ~r + —2 — ~r i A - t) i 

and, as pointed out in Sec. IV, a small imaginary part is added to the frequencies 

Ll>i = U>\r — iAj ; Ll>2 = UJ 2r ~ *-^2 ■ 

Inserting Eq. lA.5|) into Eci. HA.3fl and taking the absolute square S*S one arrives at the extremely involved result for 
-Pi2j Eq. ljA.ip . for interfering double resonances. For time reversal invariant systems which we analyze in this paper 
exclusively, one obtains 

P12M - 



£>(«) 
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with 



N(u>) 



4 { [27B,2l7B,22 A l W l,r + 2 7B 4 i7B, l 2 A 2 W2,r] + [7B ,2l7S ,22 ( A l + ^ 
+7-B,ll7S 42 ( A 2 +^ - w 2,r)] f > 




D(w) 



[ A 2(7B 2 n + 7s 2 i2 + 2AiA|wi, r ) + (7_b 2 n + 7B 2 i2) w2 + 2Aio> 2 cJi >r . 
+ (7s 2 2i + 7b 2 22)( a i + ^ 2 - ^ 2 r ) + 2A 2 (A 2 w 2 , r + w 2 w 2/r - u;i iT ,cj 2) r) 
-(7B 2 n +7s 2 i 2 ) w ir - 2AlWl >r W^ r ] 

r 2 2 r, i 2 2 \2\2/\2,\2\2 

L7B,u7B,21 ~ 2 7B,1i7B,127B,2i7B,22 + TS,ll7B,22 - A 1 A 2 _ ( A l + A 2J W 

-w 4 + 2(7 B 2 21 + 7 B 2 22 )AicJi ir + }? 2 ui{ r + u) 2 ul r + 27 B 2 n A2w2,r 

+ 2 7B 2 12 A 2W2,r + 4AiA 2 Wi, r W 2 ,r + (A 2 + OJ 2 - U^^Wj.r] 2 • 



Setting 7b 21 = 7b 22 = A 2 = w 2 , r = and dividing numerator and denominator by uj 4 one obtains the single 
resonance formula of Ea. (|31|l . 

In the parameter search of the fitting procedure one has to assure that the decay parameters Ai, A 2 stay positive. 
This can be achieved by setting Ai = (A 2 )^ 1 / 2 ^. 
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FIG. 2: Perpendicular cut through the cavity at the antenna entrances. To visualize the standing wave situation, virtual 
mirrors Mi are placed in the output antennas one quarter wavelength away from the antenna footpoints in order to produce 
the boundary conditions as expressed in Eq.l^J- I n the input antenna the mirror M a has to be placed at such a point that the 
frequency of the system stays at the continously varying value u). This reflects itself in the presence of sine and cosine waves 
in the entrance channel, Eq. 1150 , 
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FIG. 3: Comparison of two resonances measured in superconducting resonators with different quality factors. Parameters of 
upper resonance (from 0, 2D lead resonator): F 1 /2tt = 91.7 kHz, r 2 /2n = 34.9 kHz; f r = u) r /2-K = 4.63 GHz; A/2vr = 84.1 
kHz; Vtot/2-K = 294.9 kHz; 2ttt = 11.9 us; fit accuracy x 2 = 1-9. 

Parameters of lower resonance (from |T(|, 3D niobium resonator): ri/27r = 0.02 kHz, r 2 /27r = 0.18 kHz; f r = io r /2-K = 3.48 
GHz; X/2tv = 0.51 kHz; Ftot/2-iT = 1.22 kHz; 2ttt = 2.0 ms; fit accuracy y 2 = 9.7. 
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FIG. 4: Comparison of a 'warm' (T = 300 K) with a 'cold' (T = 4.2 K) cavity, showing the same isolated resonance in both 
cases. The shift of the resonance frequency is due to the reduction of the size of the cavity when cooling it to T = 4.2 K. The 
fit curves coincide completely with the measured points. 

Parameters of upper resonance: Ti/2n = 1.3 kHz, F 2 /2tv = 1.0 kHz; f r = uj r /2iv = 1.387 GHz; A/2tt = 949.2 kHz; r t ot/2-iT = 
1900.6 kHz; 2ttt = 1.1 (is; fit accuracy V = 1.07. 

Parameters of lower resonance (from Q): ri/27r = 3.8 kHz, r 2 /27r = 4.0 kHz; f T = uj r /2ir = 1.393 GHz; X/2n = 10.3 kHz; 
r tot /27r = 28.4 kHz; 2ttt = 91.1 /is; fit accuracy X 2 = 0.6. 
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FIG. 5: Two interfering double resonances (from 13). 
Parameters of upper resonance 

first resonance: Tii/2-k = 6.5 kHz, F 12 /2tv = 4.7 kHz, f lr = uj lr /2n = 2.081 GHz, Ai/2tt = 14.0 kHz, Y i>tot /2-K = 39.1 kHz; 

27TT1 = 71.4 [IS 

second resonance: r 2 i/27r = 3.0 kHz, r 22 /27r = 2.3 kHz, f 2r = aj 2 r/27r = 2.087 GHz, A 2 /2tt = 4.8 kHz, Y 2 ,tot /2tv = 14.9 kHz; 

2ttt 2 = 209.9 /is; 

fit accuracy % 2 = 0.67. 

Parameters of lower resonance 

first resonance: rn/27r = 10.5 kHz, ri 2 /27r = 9.9 kHz, f lr = Lu lr /2-iT = 10.028 GHz, Ai/2tt = 74.8 kHz, ri, tot /27r = 170.0 kHz; 
27TT1 = 13.4 us 

second resonance: r 2 i/27r = 0.0033 kHz, r 22 /2?r = 34.7 kHz, f 2r = lo2t/2-k = 10.030 GHz, A 2 /2tt = 29.7 kHz, r 2 ,t ot /27r = 94.1 
kHz; 2-7TT2 = 33.7 /is; 
fit accuracy: x 2 — 3.0. 
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FIG. 6: Part of the transmission spectrum of a 3D stadium billiard pH in the vicinity of 13.03 GHz, where the resonances 
still can be described with the fit formulae. Fit parameters for the double resonance shown in the insert, first resonance: 
rn/27T = 0.550 kHz, Fi2/27r = 0.008 kHz, f lr = Lo lr /2n = 13.0329 GHz, Ai/2vr = 1.44 kHz, r 1 , tot /2ir = 2.889 kHz; second 
resonance: T 2 i/2tt = 0.528 kHz, r 22 /27r = 0.529 kHz, f 2r = u 2r /27r = 13.0331 GHz, A 2 /2tt = 7.16 kHz, T 2 ,tot/2TY = 14.32 kHz; 
fit accuracy: \ 2 ~ 1-1- 



